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ABSTRACT 


The  energy-momentum  relationship  for  a  neutron  gas 
has  been  calculated  by  using  Hamada-Johnston ’ s  two  nucleon 
potential  and  Schwinger’s  variational  method  to  determine 
the  diagonal  elements  of  the  reaction  matrix.  An  equation 
of  state  is  constructed  so  that  the  pressure  follows 
a  simple  asymtotic  form  at  higher  densities.  It  is  known 
that  when  an  amount  of  cold  matter  is  assembled  whose 
mass  exceeds  a  certain  critical  mass  then  the  system 
becomes  unstable  against  collapse.  We  use  our  equation 
of  state  to  integrate  the  general  relativistic  equations 
of  hydrostatic  equilibrium  for  the  neutron  gas  and 
determine  a  numerical  value  for  the  critical  mass. 
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Chapter  I.  Introduction 

Recently  the  discovery  of  strong  celestial  sources 
of  X-rays  has  revived  interest  in  the  question  of  neutron 
and  baryon  stars.  In  view  of  this  it  becomes  interesting 
to  apply  nuclear  matter  methods  to  an  extended  system 
consisting  only  of  neutrons.  The  general  problem  of  a 
many  body  system  with  realistic  interactions  presents 
great  difficulties  as  far  as  exact  solutions  are  concerned 
so  by  necessity  all  techniques  will  be  approximate.  Here 
our  system  is  assumed  to  be  so  large  that  all  surface 
effects  are  negligible.  In  such  a  system  the  existence  of 
a  translational  invariance  will  simplify  the  calculations. 
We  will  assume  that  a  nonrelat ivist ic  quantum  mechanical 
description  is  valid,  that  the  nuclear  interaction  between 
two  particles  may  be  represented  by  a  static  potential, 
and  that  the  effect  of  many  body  forces  is  small.  In  the 
numerical  work  later  on  we  use  the  Hamada-Johnston  nuclear 
potential^)  which  successfully  accounts  for  the  observed 
two  nucleon  data  below  the  threshold  for  pion  production. 

Scattering  experiments  can  only  furnish  information 
about  the  forces  between  isolated  nucleons.  The  question 
remains  whether  these  forces  are  the  same  when  the 
particles  are  imbedded  in  a  nuclear  many  body  medium. 


' 


(  2  ) 

Gomez  and  co-workers  show  that  due  to  the  exclusion 
principle  the  effect  of  interaction  extends  only  to  a 
distance  small  compared  with  the  average  nucleon  spacing 
in  nuclear  matter,  and  thus  during  a  collision  two  particles 
have  only  a  very  small  probability  of  interacting  with  a 
third  particle.  The  corrections  to  the  potential 
energy  of  nuclear  matter  due  to  three  body  forces  have 

( "3  ) 

been  estimated  by  Bjorken  and  Wyman on  the  basis  of 
meson  theory  and  are  shown  to  be  about  1  Mev  per  particle. 

It  must  be  remembered  that  the  scattering  of  two 
particles  of  the  system  will  differ  from  the  scattering 
of  free  particles  because  a  particle  moving  through  the 

nuclear  many  body  medium  has  its  energy  momentum  relationship 

\ 

shifted  by  interaction  with  the  medium.  Any  theory  must 
also  take  into  account  the  fact  that  the  exclusion  principle 
restricts  intermediate  state  transitions  to  unoccupied 
states.  Thus  in  spite  of  the  fact  that  nuclear  forces  are 
strong  and  have  a  range  comparable  with  the  distance 
between  nearest  neighbours  in  a  nucleus,  the  mean  free  path 
for  nucleons  incident  on  nuclei  is  large  compared  to  the 
average  separation  of  nucleons.  The  success  of  the  shell 
model  and  optical  model  showed  that  in  some  ways  an 
independant  particle  model  is  valid.  The  reaction  matrix 
type  methods  which  we  shall  use  allow  us  to  treat  the 
particles  as  effectively  independant  provided  in  defining 
their  interactions  we  have  made  allowance  for  the  medium 
in  which  they  are  moving.  Brueckner^^  shows  that  in  first 
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approximation  nuclear  matter  is  a  low  density  system, 
so  it  is  possible  to  first  consider  the  distribution  of 
relative  momentum  by  regarding  the  system  as  an  ideal 
Fermi  gas.  This  is  a  better  approximation  for  the  neutron 
gas  as  its  density  is  half  that  of  the  corresponding 
nuclear  matter  problem.  In  particular  the  exclusion  principle 
says  that  two  particles  far  apart  will  have  insufficient 
momentum  transfer  to  raise  them  above  the  Fermi  level  so  to 
a  good  approximation  the  particles  moving  in  the  neutron 
gas  will  behave  as  independent  particles  moving  in  a  common 

C  p-  g  f-T  g  \ 

potential.  Other  authors  55  3  3  '  have  found  that  nuclear 
forces  are  not  sufficient  to  bind  the  neutron  gas  at  any 
density.  We  have  calculated  the  compressibility  of  the 
neutron  gas  and  have  found  it  not  to  differ  more  than  a 
factor  of  two  or  three  from  the  compressibility  of  an  ideal 
Fermi  gas,  so  this  approximation  seems  well  founded.  (See  Fig. 5) 
The  strong  Fourier  components  required  in  the 
description  of  the  photoelectric  effect  suggest  a  repulsive 
core  to  the  two  nucleon  potential  and  direct  measurements 
of  the  two  nucleon  correlation  function  at  high  energy  show 
that  the  effect  of  a  repulsive  core  is  clearly  present. 

This  is  usually  represented  by  making  the  potential  infinity 
positive  at  distances  less  than  about  0.5  x  10  Mm.  In  this 
case  the  matrix  elements  of  V(r)  with  respect  to  a  determinant 
of  plane  waves  are  divergent  so  ordinary  perturbation  theory 
is  not  applicable.  In  spite  of  the  hard  core  however,  we 
know  that  the  scattering  amplitudes  for  nucleon-nucleon 
scattering  remain  finite  quantities  because  the  wavefunction 


is  zero  inside  the  core.  The  radius  of  the  hard  core  is 

much  less  than  the  average  spacing  of  particles  inside  a 

-11 

nucleus  (1.7  F  where  one  F=  10  cm.).  Moreover,  the 

properties  of  the  hard  core  are  not  evident  except  at  very 
high  energies  so  we  expect  the  properties  of  the  neutron 
gas  to  depend  more  strongly  on  the  long  range  part  of  the 
potential  and  not  on  the  detailed  nature  of  the  hard  core. 

In  order  to  overcome  difficulties  associated  with  strong 
nuclear  forces,  Brueckner's  theory  replaces  the  potential 
V(r)  by  an  operator  G  which  has  finite  matrix  elements  in 
momentum  space.  The  operator  G  is  defined  by  a  non-linear 
integral  equation  which  tak-es  into  account  the  successive 
interactions  between  any  pair  of  particles  to  all  orders. 

This  integral  equation  is  similar  to  the  equation  for  the 
T-matrix  of  two  body  scattering  theory.  The  G  matrix  is 
usually  referred  to  as  the  Brueckner  reaction  matrix.  The 
equation  for  G  is  not  a  scattering  equation  however  because 
the  total  energy  which  is  defined  in  terms  of  G  itself  appears 
in  the  propagator  or  energy  denominator  and  not  just  the 
kinetic  or  free  particle  energy.  In  other  words  the  effects 
of  interaction  which  appear  in  the  denominator  change  the 
single  particle  energy  so  that  the  medium  becomes  dispersive 
with  the  energy  spectrum  (which  is  not  known  apriori)  being 
defined  in  terms  of  G  itself.  Numerical  solutions  have  been 
obtained  by  iteration^ 11 ^  but  the  requirement  of  self 
consistancy  makes  calculation  complicated  and  lengthy. 
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Some  progress  in  the  nuclear  many  body  problem  has  been 
made  using  perturbation  theories  and  rearrangement  methods. 
Although  results  seem  to  be  approximately  correct,  as  yet 
questions  on  the  existence  and  convergence  of  the  series 
are  not  clearly  resolved. (12)  Higher  order  terms  in  the 
series  expansions  have  such  complicated  structures  that 
only  the  first  few  terms  may  be  calculated.  Estimates 
of  the  errors  involved  in  any  practical  calculation  may  be 
more  tractable  when  the  problem  is  related  to  a  variational 
principle  as  we  have  done. 

It  should  be  noted  in  passing  that  the  G-matrix  becomes 
singular  for  attractive  potentials  when  the  momenta  of  two 
particles  are  nearly  opposite  and  equal  to  the  Fermi 
momentum.  The  singularity  is  a  weak  one  and  doesn’t  disturb 
computer  calculation.  The  poles  correspond  to  the  existence, 
of  bound  states  between  two  particles  and  are  connected 
with  the  phenomena  of  superconductivity.  It  is  possible 
to  ignore  this  aspect  of  the  problem  for  the  purpose  of 
calculating  energy  momentum  relations. 

The  format  of  the  thesis  is  as  follows.  Chapter  II 
contains  the  formal  developement  of  the  equations  to  be 
solved  for  the  reaction  matrix  and  related  quantities  such 
as  energy  per  particle.  Chapter  III  relates  the  problem 
to  a  variational  principle  and  contains  the  methods  and 
results  of  numerical  solution  along  with  a  comparison 
with  other  authors.  Chapter  IV  is  an  introduction  to  the 
subject  of  superdense  stellar  configurations  and  a 
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derrivation  of  the  relativistic  equations  of  hydrostatic 
equilibrium  from  Einstein's  field  equations.  Chapter  V 
combines  the  results  of  chapters  III  and  IV  by  integrating 
the  equilibrium  equations  using  the  equation  of  state  we 
have  developed  for  the  neutron  gas.  We  discuss  the  question 
of  stability  of  the  star  and  give  a  numerical  value  for  the 
critical  mass  beyond  which  there  are  no  stable  configurations 
of  a  neutron  star.  The  question  of  what  happens  to  a  star 
at  the  endpoint  of  thermonuclear  evolution  which  has  a 
mass  greater  than  the  critical  mass  is  discussed  briefly. 


7. 


Chapter  II.  The  Nuclear  Reaction  Matrix 

The  total  Hamiltonian  for  a  system  for  N  non-relativist ic 
particles  may  be  written; 

H  =  K-tfv?)  +  1  l  V(r.,r,)  (1.1) 

1  2m  2  JYj  J 

ie .  the  sum  of  kinetic  energy  contributions  and  an 
interaction  part  involving  only  two  body  interactions  (we 
are  assuming  that  three  body  forces  may  be  neglected) . 

Write 


H  =  H0+  H'  (1.2) 

where  H'  represents  the  total  interaction  term.  Ho  defines 
a  set  of  wave  vectors  |  $j_>  by 

H0 |$±>  =  Ei|$i>  (1.3) 

If  there  are  N  particles  in  the  system,  then  $  is  a 
determinental  wave  function; 

$  =  (N!  )-%  Det  |  (J>  i  <}>2  •  •  •  (|>N  |  (1.4) 

where  cf^are  single  particle  wave  functions  defined  by 

H  o  4>  =  e±$±  (1.5) 

$o  is  the  state  obtained  by  putting  the  particles  in  the  N 
lowest  available  levels  <J>j_  .  H'  represents  the  two  body 
interaction  between  particles  which  may  either  change  the 


' 


■ 


state  of  two  particles  or  none.  Singly  excited  states 
cannot  be  reached  because  of  momentum  conservation.  In 
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the  representation  defined  by  the  set  {$},  the  two  particle 
matrix  elements  of  an  operator  V  are; 


<r  s  |  V  |  mn>  =  <fr  ( 1 )  cj>s*(  2  )  Vi  2<J>m( 1 )  <f>n(  2  ) dx  i dx  2  (1.6) 


For  an  infinite  system  the  basic  single  particle  wave 
functions  are  plane  waves. 


(1.7) 


where  ft  =  L3  is  the  large  volume  containing  the  N  particles. 
Periodic  boundary  conditions  imply  that  each  component  of 
It  can  take  on  values  27m/L;  (n  =  0,  ±1,  ±2,  •••)•  In  the 
limit  that  ft->-°°  while  N/ft  remains  finite  and  constant, 
summations  must  be  replaced  by  integrations.  i.e . 


(1.8) 


(there  are  states  per  unit  volume).  In  the  ground 


state  of  the  non-interacting  system,  the  particles  occupy 
a  sphereof  radius  kp  in  momentum  space  and 


kp  -  p  ( kp ) 

3  TT  2 


N 

ft 


(1.9) 


where  p(kp)  is  the  particle  density. 

The  total  Hamiltonian  H  specifies  the  wave  vector  | I 0 > 
and  energy  E  of  the  system  in  the  ground  state. 


H  ?0> 


E  hF  0  > 
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Expand  \'¥q>  in  terms  of  the  complete  set  |$n>,  ie . 


fo>  =  l 
m 


|  T0  > 

=  l  anl*n 

> 

n 

(H 

o  +  H')I  an 

brC  =  EI 

a  1 
n  1 

0  > 
n 

n 

n 

I 

anEn | $n>  + 

1  anE  ’  | 

= 

El  a 

n 

n 

n 

amEm  +  1  an<^,m|H,|<^n 

>  = 

Eam 

n 

am  -  y 

an<(^m  |  H  1  |  0 

n> 

n 

E-Em 

m  1 

*m>  =  l  hm 

>  jl  < 1  H 1 

1  ^n 

4 

> 

+ 

m^O 

E  - 

Em 

(1.10) 


n' 


+  aoN>o>  (1.11) 


The  state  |O0>  is  normalized  to  unity  but  |  T  0  >  is  normalized 
so  that  its  scaler  product  with  |<f>o>  is  unity.  Then, 


y„>  =  |4>0>  +  l  I  <t,m>Q<  “I'm  I  H '  po? 

m^O  - 


(1.12) 


E  -  E 


m 


also 


E  =  E0  +  <$o  |  H  '  I  T0  >. 


(1.13) 


A  projection  operator  Q  which  excludes  occupied  states  has 
been  inserted  to  take  care  of  the  exclusion  principle. 
Assuming  that  the  center  of  mass  momentum  of  two 
interacting  particles  is  zero,  Q  will  be  unity  if  in  the 
intermediate  states  both  particles  are  outside  the  Fermi 
sea  and  zero  if  one  or  both  are  inside.  Q  depends  on  the 
relative  and  total  momentum  in  the  intermediate  states  so 


for  non-zero  center  of  mass  momentum  it  will  have  a  slightly 
different  form.  Operate  on  equation  (1.12)  by  e  =  E  -  Em 
and  note  that  since  |$0>  is  a  plane  wave  state  e|$0>  =  0. 

In  operator  form  (1.12)  becomes 


el o  =  QVI o 


(1.14) 


which  demonstrates  its  similarity  with  a  two  particle 
Schrodinger  equation  (with  a  non-local  potential). 

For  each  pair  of  particles  define  a  wave  vector  1 1  ,•  •  > 
which  contains  correlations  of  all  pairs  except  (i,j)  by 


y  >  =  | $  0  >  +  l  | 0n>Q< $n | H ’  -  V •  •  1 1 o  > 

n^O - s - 


(1.15) 


E  -  E 


n 


Define  a  two  body  operator  G •  •  by 

J 


Gij I ^i  j >  vij I >  ' 


(1.16) 


Equations  (1.12)  and  (1.15)  combine  to  give 


Io> 


+  l  |  ■fn>Q<tn|Vlj.  |T„> 


(1.17) 


Operate  on  (1.17)  by  Vj_j  and  assume  |^ij>  -  |$o>  which  is 
equivalent  to  neglecting  many  particle  excitations.  Then 


<  $m I Gi j  I  $ o  > 


<$m|V 


ij 


$o  > 


+  I  I  ^i  i  I  (^n>Q<^n  I  Gi  i  I  $  o  > 

n^O  - ^ - ±4 - .(1.18) 

E  -  En 


Io> 


We  have 


E  =  E0  +  <$o|  H  ’ 


11. 


or>  E  =  E0  +  h  l  <$o  |V.  ,  \V0> 

ij  J 

=  E  o  +  3g  l  .<$0  |  G  j_  j  \y±j> 


~  E0  +  h  l  <$o  | G  | $o>  (1.19) 

ij  1J 


Since  the  matrix  elements  of  G-,  depend  only  on  the  • 
single  particle  wave  functions  of  particles  i  and  j ,  the 
integral  equation  for  G  becomes, 


<i '  j  '  I G  |  ij >  =  <i'j'  |  V | 


J 


mn^ij 


ei+  ej -  em-en 


(1.20) 


where  ej_  represent  single  particle  energies  defined  in 
terms  of  G  itself.  These  will  be  discussed  in  detail  later. 

It  is  conveniant  to  introduce  the  so  called  Mdller 
wave  matrix  ft  by 


V  •  •  ft  •  • 


(1.21) 


Equation  (1.20)  in  the  coordinate  representation  for  two 
particles  in  states  and  Z2  becomes 


->  -> 


<r i r 2 | ftp j | £ i £2  >  =  <rir2|£i£2>  +1  <r i r 2 | £ f &2 >Q<£ f £2 ' V± . ft . , | £ x &2> 


£|  1'2 


iJ  ij 


e£l  +  e  1 2 


"  e£;  “  et 2 


let 


<r x? 2 | ft | £ i £2  >  E  Y£l£2(?i ,?2) 
<r ir2 | ft | £i £2>  E  $£l£2(?i,r2) 


(1.23) 


* 

* 


* 
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G(?i  ,r2  ,r  2  )  =  l  £  -  ( ?i  ,r  2  )  Q$l  f  £ f  ( r  T  ^2  ) 

£  1  2/  2 - 


* 


£1 


+e 


£2  ~e£r  ~e£ 


<?i srI I  V | r 1 ,r 2>  =  V(?i,rT ;?!,?$) 


(1.24) 


Then , 


-*  -> 


J.jJ.fruU)  =  4'£1£2(?1?2)  + 


G(r i?2 ;r r?2 )V(r r?2 


X  £2  (rTr~)dr  Tdr^dr  "dr" 


(1.25) 


If  V  is  local  and  doesn't  depend  on  the  position  of  the 
center  of  mass  then, 


<r ir2 | V | r r?2>  =  6 (r 1 +r2-r 7-r 2 ) 6 (r 1 -r 2-r f+r 2 )V( | r 1 -r2 | ) 


-y 


(1.26) 


-y 


llSi2(r1rz)  =  £2  (r  >r2 )  + 


G(r  ir2  ;rlrl  )VT  (rfr£)drfdrj 

X/  I  Xj  o 

(1.27) 


We  are  considering  a  system  of  infinite  extent  so 
$(rir2)  are  plane  waves  and  therefore  separable  in  terms  of 
center  of  mass  and  relative  co-ordinates  R  =  ^(rj  +  r2) 
and  r  =  ri  -r2.  Write  the  quantum  numbers  for  center  of 
mass  and  relative  motion  corresponding  to  £1  and  £2as 
k  and  P .  Then 


■y  -y 


-y  -y . 


^kP ( r  s R )  =  ^kp(^,R)  + 


kP 


G(r  ,R;r;RjV(r")dr  "dR"  Y  (»',(?)  (1.28) 

1  ■ 


The  $'s  are  separable,  e.g . 


*kp(?,S)  =  l  4k(r)*p(S) 


(1.29) 


i 
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Then,  ¥kp(r) 


(1.30) 


and  Gp(r,r’)  =  l  $k , ( r ) $k » (r ’ ) Q( ? ) 

k  ’  - 

ekP  ~ek’P 


(1.3D 


Equation  (1.28)  is  the  standard  form  of  the  Bethe-Goldstone 
equation.  The  matrix  element  of  G  is 


<£’ I G |£> 


dr  ’  $k»  (r  ’  )V(r  ’  )'tkP(r  ’  ). 


(1.32) 


As  is  usual  in  the  two  body  scattering  problem,  we  attempt 
to  make  a  partial  wave  expansion  following  the  procedure  of 
Brueckner  and  Gammel.^^  It  is  conveniant  at  this  point  to 
remember  that  we  must  consider  the  particle  spin  for  a 
complete  specification  of  the  dynamical  state  of  the  system. 
A  complication  arises  here  because  in  equation  (1.31)  for 
the  Green’s  function  the  summation  depends  on  the  angle 
between  ?  and  k’  through  both  the  energies  appearing  in  the 
denominator  and  through  the  exclusion  operator  Q(k’,P) 


Q(k 


=  0 

for 

|  k’ 

i+ 

l_cM' 

£  kp 

=  1 

for 

|£» 

+  ?| 

>kp 

(1.33) 

To  remove  this  coupling  it  is  usual  to  perform  an  average 
over  angles.  In  particular  the  averaging  of  Q  is  a  good 
approximation  as  the  G-matrix  seems  to  depend  only  weakly 
on  the  form  of  the  exclusion  principle. 


14. 


Q(k',P)  =  0  ( k '  2+P2+2k ’  Pa-k|i )  0  ( k  ’  2+P2-2k ’  Pa-kp )  (1.34) 


where  a  =  cos(ic’,?) 
Put  a  =  k» 2+P2-kp , 

Define  Q(k’ ,P)  =  \ 


,  and  the  0's  are 

and  b  =  k  '  P 

+  i 

* 

Q(£' ,?)da  =  % 


step  functions. 

+  i 

0 ( a-ba ) 0 ( a+ba )da 


(  .35) 


y  - 1 


'  -  i 


Q(k* ,P)  =  0 ( a-b ) 0 ( a+b )  +  a  0(2a)  (1.36) 

b 


ie .  Q  =  0  for  k’ 2+P2  <kp 

=  1  if  |k’!-|?|  >kF  (1.37) 

=  k,2+P2-kp  otherwise 

2k  ’  P 


The  angular  dependance  in  the  denominator  may  be  seen  by 
expanding  the  single  particle  energies  for  exicited  states 
in  powers  of  momentum. 


em  +  en  =  2Uo  +  Ui-(k^  +  k£)  +  U2'(k^  +  k^)  +  ...  (1.38) 

-  2U o  +2Ui*(k’  2+P2  )  +2U2'(k,4+2(ic’  *?)2+P4+k,2P2)  +... 

The  angle  between  k'  and  ?  occurs  first  in  the  quartic 
term.  This  could  be  eliminated  by  taking  an  angle  average 
of  the  term  (k'*?)2  .  Bethe^1^)  has  argued  to  neglect 
the  single  particle  potentials  in  the  excited  states  for 
a  pure  neutron  gas.  We  follow  his  example  and  proceed  with 
the  partial  wave  expansion.  The  non-interacting  two 
particle  wave  fucntion  for  relative  motion  may  be  expanded  as 


* 


where 


f 


=  ( 2£+l ) i 


£ 


4 


2  £+1 


and  f(R,f)  means  a  function  of  the  angles  between  k  and  r. 
In  the  presence  of  a  tensor  force,  the  orbital  angular 
momentum  number  £  is  not  a  constant  of  motion  but  the  total 
angular  momentum  is.  Equation  (1.39)  should  be  expanded 
in  terms  of  eigenfunctions  of  J,  that  is  Fl™! S .  These 
functions  are  related  to  the  above  by. 


pJmjs 

L 


+i 


=  I  Yl  ms  Xss  C(J,mJ;L,mJ-ms;s,ms) 


(1.40) 


ms  =  "1 


The  C’s  in  the  above  expression  are  Clebsch-Gordan 
coefficients.  Inverting  equation  (1.40), 


LV1UJ  jmL+ms  > s 

l  FL 

J  =  L-1 


C( J ,mL+ms ;L,mL; smg ) . 


(1.41) 


Equation  (1.39)  then  becomes, 

00  J  +  l  jm  s 

tk(r)  =  l  l  c£  j^(kr)  P£  s  (R,f)  C  ( Jms ;  £0 ;  sms )  (1.42) 

J=0  £=J-1 


We  want  a  similar  expansion  for  T^(r).  It  must  be 
remembered  here  that  the  tensor  force  couples  the  states 


for  s  =  1,  £  =  J  ±  1,  (the  state  s  =  1,  £  =  J  is  not 
coupled  with  any  other).  We  adopt  a  notation  for 
distinguishing  the  states  £  =  J  ±  1  by  putting  subscripts 
on  the  partial  wave  functions  U. 


No  tensor  force 


With  tensor  force 


J±1 


,  s  Jm j  s 
(kr)  PJ_1  (k,f) 


U  (kr)FJmjS+U  JmjS 

UJ±1,J-1U  ^J-l  +UJ±1,J+1FJ+1 


(1.43) 


Denote  the  uncoupled  triplet  states  by 

JniTS  Js  JmTs 

jj(kr)  Fj  -  Ujj ( kr )  Fj 


Then  the  expansion  for 

00  J  +  l  J  +1  T  o 

v?>  -i  i  i  At 

J=0  £  J_1 £ '=J -1 


T^(r)  is  (dropping  the  index  P  on  V ) , 
(r)  F^sS(k,r)  C  ( Jms  ;  £0  ;  sms  )  .  (1.44) 


A  similar  expansion  for  the  Green’s  function  is. 


G(r,r")  =  l  c£i^  G^(r,r")  Yj(r,r")  (1.45) 

£  =  0 


where 


G£(r 9v') 


1 

2  7T  2 


k"2dk"j£(k"r)  j£(k"r")  Q(k >P) # (1 e 45 ) 


0 


e(k,k;p) 
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Insert  these  expressions  into  equation  (1.28)  and  drop  the 
sum  over  J  which  is  common  to  all  terms. 


„  Js  Jmss 

l  uJU'(r)  (&»£)  C(Jms;£0;sms) 


Jmss 


I  C  j  £ ( kr )  PA”B  (K,f)  C(Jms;£0;sms) 


£ 


+  I 
£  > 


dr "c£(-i) £G£(r 3r')  Yjtf^f')  V(r") 


x  l  c 


Z'Z" 


Js  JmQs 

t'Un^.lr')  F£„s  (R.f)  C(Jms;£0;sms) 


(1.47) 


Integrate  over  the  angles  of  r’ 


/  .  \  Z 

(~i)  c 


Z 


n  Jrrio  s 

dr'  Y°(f>  ,f"  )  V(r '  )  F£„s  (R,f>0 


Jm^s  Js 

4 it  F£  S  (R,f»)  Vu„(r') 


(1.48) 


where  V^£M(r)  = 


Jm0s  Jm.-.s 

df>  F£  s  (R,r)  V(r)  F£„s  (R,f) 


(1.48) 


and  we  have  used  the  relationship 


Y ° (r ,f  ’  ) 


r  4 


2  £+1 


i 


l  Y£(R,f')  Y"(R,£>) 


m=-l 


£ 


(1.50) 


Js  Js 

V££"(r)  =  spp"  vT^r^  for  tensor  forces 


££ 


=  Soon  v„(r)  for  central  forces 
£  £  ^ 


We  may  write 


18. 


Here  vT(r)  and  vc(r)  are  radial  dependances  and  the  S’s 

are  numbers  listed  by  Ashkin  and  Wu^^),  Put  equation 

/  i  \  /I\  JmQs 

(1.48)  into  (1.47),  take  the  scaler  product  with  F  f s  (K,r), 

integrate  over  angles  of  r  (change  the  dummy  index  £  to  1' 

on  the  left  hand  side).  Then  we  get, 


I 

£ 


Js 


(r ) 


C ( Jmg ; £ ’ 0 ; sms ) 


CL , j  L , ( kr )  C ( Jms ; L ' 0 ; sms ) 


+  4tt^  J  ^^r  ’  2dr  '  GL,(r,r')  V^£„(r)  c0,  0  „  (r  '  )  C  ( Jms  ;  £0  ;  sms  ) 


'£ 


£  »  £' 


(1.51) 


Now  multiply  through  by  C(Jms;L0;sm  )  and  sum  over  m 

^  o  o 

using  the  orthogonality  relationship 


1 

l  C(Jms;L0;sms)  C ( Jmg ; £ ' 0 ; smg ) 


L£ 


2  J  +  l 
2  £  ’  +1 


(1.52) 


which  gives  (changing  L+£  and  L’-H’) 


Js 

Uoo.  =  JVkr)6££  +  4tt  ^ 


r ' 2dr’G0 , (r,r! )V0  0„(r ’ )U0  0„(r T  ) 


£ 


££ 


££ 


£" 


(1.53) 


The  equations  for  all  angular  momentum  states  have 
similar  form  so  let  us  suppress  indices  for  a  moment. 


U  ( r ) 


U  o  ( r )  +  4tt 


G(r,r’)  V(r’)  U(rf)  r,2dr’ 


(1.54) 


- 
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This  must  be  modified  for  potentials  with  hard  cores  because 
V(r)  =  oo  over  the  core  region.  Brueckner  and  Gammel*'11^ 
show  for  the  two  body  scattering  problem  the  replacement 


V( r )U(r )  =  A6 ( r-c ) 


(1.55) 


for  r  ^  c  leads  to  the  correct  Schrodinger  equation  for  two 
body  scattering.  A  is  determined  by  the  boundary  condition 
U(c)  =  0. 


U(r)  =  U0(r)  +47Tc2AG(r  ,c )  +^tt 


G(r ,r  1 )U(r ’ )V(r f )r'  2dr ' 

(1.56) 


U(c)=0  implies  A= 


G(cJr,)V(r')U(r')r,2dr’  +U0(r)  (1.57) 


G  (  c  ,  c  ) 


Equation  (1.5*0  is  replaced  by 


U(r )  =  f  (r )  +  4tt 


K( r ,r  ' )  V( r ' )U(r ’ )r,2drT  (1.58) 


with  K(r,r’)  =  G(r,rT)  -  — —  ^-c  }  r — — 

G  (  c  ,  c  ) 


and  f  (r )  =  U0(r)  -  U0(c)  G-(-c-^-r>-) 

G  (  c  ,  c  ) 


For  the  many  body  system,  the  situation  is  not  as  simple. 

( 22 ) 

However,  Bethe  and  Goldstone  show  that  replacing  equation 
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(1.5*0  by  (1.58)  and  setting  U(r)  =  0  for  r<  c  leads  to 
negligible  error.  Brueckner  and  Gammel^11^  have  also 
considered  this  problem  for  a  hard  core  potential  with  no 
attraction,  taking  S-states  only.  They  write 

V(r)U(r)  =  l<$(r-c)  +  f(r)  for  r  <  c 


and  then  show  the  effects  of  the  correction  term  are  about 
0.1%  for  nuclear  matter  at  equilibrium  density.  We  therefore 
feel  justified  in  following  the  same  approach  for  a  neutron 
system.  The  general  treatment  taking  into  account  angular 
momentum  states  is  exactly  the  same  as  was  just  done  and 
consists  of  making  the  following  replacements  in  Equation 
(1.53;  . 

Gf(c,r) 

j£(kr)  +  f^(kr)  =  j^(kr)  -  j£(kc)  _  (1.60) 

G^(  c  ,  c  ) 

G^(r,r'  )  K^(r,r’)  e  G^(r,r’)  -  ?c  ^  c  ?r  ^  (1.6l) 

G£( c ,c ) 


Inserting  equations  (1.42)  and  (1.44)  into  (1.32),  we  find 
for  the  diagonal  elements  of  the  reaction  matrix, 


J+l 

=  4tt£  l  (  2  J+l ) 
J  £=J -1 


dr j  ^ ( kr ) 
£ 


7*1 


l 

'=j-i 


yJs  u^S 
V  £#  ’  U££  ’ 


<k  G  k> 


OO 


(1.62) 
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This  expression  may  be  separated  into  two  contributions; 
one  representing  the  core  volume  term,  the  other  the  effect 
of  the  core  surface  and  outer  potential.  The  procedure 
is  the  same  as  before.  Replace  V(r')U(r’)  by  A6(r’-c) 
for  r'  ^  c,  where  A  is  given  by  (1.57)  and  the  core  radius 
is  assumed  the  same  in  all  states.  Then, 


<k  G  k>  =  <k  G  k> 


out 


+  <k  G  k> 


core 


oo 


J+l 


+  r 2dr  f£(kr)  J  V^®(r)  U 


£ ’ = J-l 


^  (r)  (1.63) 


+  <k  G  k> 


core  ’ 


The  contribution  from  the  core  was  calculated  from  equation 
(2,3)  of  reference  16. 


. 

- 


*,  ^  K>, 
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Chapter  III.  The  Equation  of  State  for  a  Neutron  Gas 

In  attempting  a  particular  solution  of  the  problem  for 
an  extended  system  consisting  only  of  neutrons,  we  must 
remember  that  the  exclusion  principle  states  that  the 
wave  function  must  be  antisymetric  under  exchange  of  all 
particle  co-ordinates.  For  the  neutron  gas  we  will  have 
only  to  consider  singlet  states  of  even  parity  and  triplet 
states  of  odd  parity.  We  believe  that  the  system  will  be 
adequately  described  by  considering  only  the  partial  waves 
for  i  =  0  and  1.  i,e .  the  contibutions  to  the  reaction 
matrix  for  i  >  1  are  very  small. ’  We  make  the  further 
assumption  here  that  the  coupling  between  the  3P2  and  3F2 
states  is  sufficiently  weak  that  we  may  consider  the  tensor 
force  for  the  3F2  state  to  be  identically  zero. 

Equation  (1.53)  can  be  solved  only  by  making  certain 
approximations  whose  validity  is  difficult  to  assess. 

We  now  attempt  an  approximate  solution  by  using  a 
variational  technique  similar  to  Schwinger  ’  s  .  ^^  ) 

Consider  equation  (1.53)  in  the  form 


U£>  =  lf£>  +  K£(r,r')  V(r’)|U£>.  (2.2) 


The  reaction  amplitude  for  the  outer  potential  is 


%  ■  <fJvlV- 


(2.3) 
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From  equation  (2.2)  we  may  write 


<uJyIV  -  <uJvVlV 


<fJvlV 


or 


<fJvIV  “  l<fJyIVl2 


<UJVIU£>  -  <U£lVK£VlU£> 


(2. it) 


Consider  the  quantity 


J(Unm)  =  I <f  „ I V| U„m> 


£T 


i 


IT' 


(2.5) 


<U£T I V I U£T>  “  <U£T 


VK„V  U.m> 


■£ 


AT' 


where  | U£T>  is  some  trial  wave  function  satisfying  the 
same  boundary  conditions  as  |U^>,  namely  U(^r+°°)  j^(kr) 
and  Ulc^  =  ^  where  0  ls  the  radius  of  the  hard  core. 

The  quantity  J(U^T)  has  the  property  of  being  homogeneous 
I U£T>  so  if  we  multiply  |U^>  by  some  arbitrary  constant 
factor,  the  expression  remains  unchanged.  Note  also  that 
lu£T>  is  important  only  if  it  is  within  the  range  of  nuclear 
forces.  Also,  the  error  in  J(U^^)  is  proportional  to  the 
square  of  the  error  in  |Uom>.  Write 


J  ( U p n-i )  =  P2  .  where 

Q  5 


P  =  <f£|V|U£T>  (2.6) 

Q  -  <UAT|V|UJtT>  -  <U£T|VKV|UiT> 


2b . 


Consider  the  variation  of  J(U^T)  with  respect  to 
variations  in  |U^^>  subject  to  the  condition  that 
satisfy  the  same  boundary  conditions  as  |U^>. 


arbitrary 


U£T> 


J  =  2P_.  6P  -  P^6Q  where 
Q  Q2 


6p  =  <f£|V|6U£T> 

6Q  =  2<U£T|V-VKV|6U£T> 


6  J  =  2P 

Q"2" 


<f „  V 


£ 


6U£T><U£T 


V-VKV 


U£T> 


-  <f  J  V 


£ 


U£T><U£T 


V-VKV 


<SU£T> 


2P 

Q2 


<U£T | V-VKV | U£T><ft | V-<f £ | V | UAT><U£T | V-VKV 


> 


=  2P 
Q2 


-  {<v 


fiU£T>} 


(2.7) 


where  <F„  =  <f„  V  U„m> 


£ 


£ 


£T 


<U£T 'V-VKV|Unm> 


£T_  .  <f  ^  |  V-<U0T,  |^7-VKV3 


£T 


<f  „  V  U„m> 


£ 


£T’ 


(2.8) 


For  6J  =  0  it  is  necessary  that  <F^  |  =  0,  i.e .  necessary 
that  |U^>  satisfy  the  integral  equation  (2.2).  Note 
that  if  lU£T>  =  |U^>,  then 

J(U£T)  =  <f£|V|U£>  (2-9) 


which  is  the  contribution  of  the  outer  potential  to  the 
reaction  matrix.  It  is  hoped  that  if  the  trial  wave 
function  is  well  choosen,  then  varying  some  arbitrary 
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parameter(s)  to  find  a  stationary  value  of  J(U^T)  will  give 
a  wave  function  that  reasonably  resembles  the  ground  state 
wave  function  of  the  system.  We  have  choosen  a  variational 
approach  to  this  problem  because  it  will  be  possible  to  put 
a  bound  on  the  ground  state  energy  of  the  system  according 
to  the  general  statements  that  can  be  made  of  variational 
principles.  The  following  considers  this  problem.  Ja  is 
the  quadratic  functional  defined  by  a  potential  Va  and 
arbitrary  | U^T> , ( sat isf ying  the  boundary  conditions  of 
course ) . 


ja  = 


<f£lvalU)lT>l2 


<UAT|Va|UJlT>  -  <U£T|VaKVa|UjlT> 


(2.10) 


Consider  another  potential  =  Va  +  e  (2.11) 

where  e  is  infinitesimal,  real,  and  positive.  Define  Jb  by 

Jb  =  _ I  <f£  lVb  I  2 _  (2.12) 

<U£T I Vb I U£T>  "  <U£T I VbKVb I U£T> 

Jb  =  I  <f£  I  Va  I  Ui^T>  I  2+  <f  ^a  lU£T><U£T^f£>  +°(£2)  • 

<U£T I Va I U£T>_<U£T I VKV I U£T>+2£<U£T I U£T>_2eU£T I VaK I U£T> 

Write  =  A  ,  and  =  A  +  eC  (2.13) 

a  B  B  +  eD 


where  A  =  I <f £ I va I U£T> I  2 

B  ■  -  <U£TlVaKValhT> 

C  =  2<f4|Ta|UtT><UJlT|fJl  > 

D  “  -  2<UATlVaKlD£T> 

Expand  Jfe  =  Jb(e=0)  +  eJb’(e=0)  +  0(e2)  +  ... 

Jb,(e=0)  =  BC  -  DA  ;  Jb(e=0)  =  J& 

B 


Therefore 


(jb  -  Ja) 


BC  -  DA 


Assume  that  Va  is  an  attractive  potential,  ie. 

hJhh  ■  -W°lhtT>  5  W°  £  0 

Consider  BC  -  DA, 

=  (-W„<U£T|U£T>-W?<UW|K|U(,T>}  (-2W0)  |<f£|U£T> 
-Wo | <f£|U£T> | 2{<u£T!u£T>  +2W0<U£t|K|U£t>} 

W°<U£T  I  U£T>  I  <fi  I  U«.T>  C 

■  =  0 


(2.13) 


(2.14) 


(2.15) 


Therefore,  Jb-Ja 


(2.16) 
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The  requirement  that  e  be  infinitesimal  can  be  dropped  as 
any  finite  difference  in  potential  may  be  constructed  from 
a  number  of  infinitesimal  ones  with  the  corresponding  shifts 
in  J.  We  may  now  write  the  following  theorem. 

1 .  For  attractive  potentials 

If  <4>  |  |  <j)>  1  <c|)|Va|(t>>  for  all  <j> ,  then  =  Ja. 

2 .  For  repulsive  potentials 

If  <(j)  |  |  (|)>  £  <cj)|Va|(f)>  for  all  <J) ,  then  Jb  =  J&. 

It  will  be  conveniant  now  to  consider  a  different 
normalization  of  the  wave  function.  Define 


K ’ (r ,r ’ )  =  K(r,r’)  V(r’),  and  Q  =  V(1  -  K?) 


where  is  a  real  Hermitian  operator.  Then 


J 


<f„  V  U„m> 


Z 


ZT' 


<U 


ZT 


Q  U 


£T 


> 


Assume  <c f  ^  |  V  |  U > 


1,  then 


1 


<U 


ZT 


S2 


UZT> 


J(uZT) 


(2.17) 
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For  an  attractive  potential  V(r),  Schwartz’s  Inequality 
states  for  any  <p  and  U^T; 


<+|V|(t»  4  <4.|y|U£T><U),T|vU>  =  <  <t>  I  Vs  1 4>  >  » 

<U£T  I V  I  UJtT> 


(2.18) 


Where  Vs  is  a  non-local  operator.  Consider  |irs>  the  wave 
function  defined  by  the  non-local  potential  and 


|UVs>  =  |f  >  +  KV  |Uvs>. 


V, 


£ 


(2.19) 


Then  <UpJV|UVs>  =  <U£,rp|  V|  f ,  >  +  <U_  |  VKVS  |  Uv  s>  ,  or 


V< 


£T 


£T  1  1  £ 


£T 


<U0T|V|UVs>  =  1  +  <Uorp|VKV|Uorp><Uorp|V|Uvs> 


v£ 


£T 


£T 


£T  £T 


<U£TIVIU£T> 


£T 


£T' 


<U£T I V I U£T>  "  <U£TIVK,IU£T>  =  <U£TIV1U£T> 

<U£T I V I UVS> 


=  <U£TI^IU£T>  E  — 1 


J(U£T} 


(2.20) 


The  reaction  amplitude  for  the  non-local  potential  Vc 
is  given  by 


<f£|Vs|UVs>  a  <ft|v|UtT><UtT|V|UVs>  =  J(UtT)  from  2.20. 


£  1  1  £T  £T 


£T 


<U£tIVIU£T> 


- 

. 

r 

■ 

>• 


29. 


Any  |U^T>  that  we  choose  will  be  equivalent  to  picking  some 
Vs  satisfying  (2.18).  J(U£T)  will  be  the  reaction  amplitude 
for  the  potential  V  ,  and  from  the  previous  theorem  (for 

o 

attractive  potentials). 


J(U£)  £  J(U£T)  (2.21) 

For  an  attractive  potential  V,  the  quantity  J(U£T)  provides 
an  upper  bound  for  the  reaction  matrix  due  to  the  potential 
V.  Similarly  for  a  repulsive  potential  the  quantity  J(Uarp) 
provides  a  lower  bound. 

For  any  practical  solution  we  are  still  faced  with  the 
fact  that  the  energy  denominator  appearing  in  the  Green’s 
function  G^(r,r')  depends  on  the  reaction  matrix  itself. 

It  contains  both  the  kinetic  and  potential  energies  of  the 
initial  and  intermediate  states.  (see  equation  1.20) 

e(k,k',?)  =  e-j_  +  ej  -  em  -  en 

=  e(  —  +£)  +e(~  -k)  -e(  —  +k '  )  -e(~  -k’) 


Here  e^  are  single  particle  energies, 

ei  =  ki  +  U(k. ) 

2m  1 


and  U(k^)  are  single  particle  potential  energies  seen  by 
a  particle  of  momentum  k^ . 


30. 


(14') 

For  the  pure  neutron  gas,  Bethev  '  has  argued  to  neglect 


the  single  particle  potential  energy  in  the  excited  states. 
Then 


For  particles  inside  the  Fermi  sea 


00 


(2.24) 


where  N(k_,k)dk  is  the  number  of  particles  in  a  degenerate 
Fermi  gas  of  N  particles  which  have  a  relative  momentum 
between  k  and  k+dk  with  respect  to  a  given  particle  of 
fixed  momentum  k^ .  This  function  may  be  obtained  as  follows. 


FERMI 


SPHERE 


For  a  given  k^,  the  maximum  of  k  is  3g(k+k^).  The  factor 
of  two  appears  because  the  relative  momentum  is  yv  where  v 
is  the  relative  velocity  and  y  is  the  reduced  mass  =  m/2. 
For  k  <  3g(kp-k^)a  k  can  go  over  all  directions.  Then 
N(kiSk)dk  is  the  number  of  particles  in  a  spherical  shell 
with  radius  2k  and  thickness  d(2k). 


31. 


For  this  case, 

N(k_  ,k)dk  =  N  4tt  (  2k)  2d  (  2k )  =  N  24k2  dk.  (2.25) 

(  4  tt/  3  )  kp  kp 

For  k  >  ^(kp-kp),  0  cannot  go  into  the  solid  angle  2tt(1-cos0) 
From  the  diagram. 


cos ( 0 ) 


4k  •  k . 


i 


Then , 


N(k  ,k)dk  =  N(2k) 2d(2k) 

(  4tt/3  )  kp 


4  7T  —  2  7T  (  1-COS0  ) 


=  N • 12 ( l+cos0 ) 


(2.26) 


The  total  energy  of  the  ground  state  of  the  system  may 
be  expressed  in  terms  of  the  diagonal  elements  of  the 
reaction  matrix  by  taking  matrix  elements  of  the  total 
Hamiltonian  (1.1)  between  the  states  |$0>  and  |4,0>  with 
<$o |¥o>  =  1. 


E  =  <  $  o  |  H  |  ¥  o  >  =  l  <  $  o  |  T  .  |  ¥  o  >  +  hi  <$o  |V.  .  |  'F  o  > 

i  ij  1J 


=  l  +  h  l  <ij | G | ij  > 
i  2  ij 


(2.27) 


32. 


The  average  of  the  kinetic  energy  term  over  all  occupied 

—  3  h-2 

levels  of  the  degenerate  Fermi  gas  is  T  =  ^  . 

The  average  potential  energy  of  the  particles  in  the  ground 

state  may  be  written  in  terms  of  the  single  particle 

potentials  or  by  finding  the  average  value  of  G  when  kp  and 

k.  take  on  all  values  in  the  Fermi  sea. 

J 


U  = 


U(k. )  d  3k . 
1  i 


k .  <kD 
i  F 


d  k . 


k .  <k„ 
i  F 


1 


P(k,kp)  < k | G | k >  dk  (2.28) 


Here  P(k,kp)dk  is  the  number  of  pairs  in  a  degenerate  Fermi 
gas  of  N  identical  particles  with  relative  momentum  between 
k  and  k+dk.  This  function  may  be  obtained  from  N(k^,k)dk 
by  integrating  k^  over  all  possible  values  (and  normalizing). 


P(k,kp)dk  = 


kF 


» 

N(k.  , k )dk •  ^Trk?dk . 

i 5  li 

=  24Nk2dk 

i  3 

fl-3k 

2k„ 

+-(-  )3} 
2^k  '  { 

(4tt/3)  k 3 

kF 

F 

F  j 

For  low  enough  densities  it  is  known  that  the  diagonal 

elements  of  the  reaction  matrix  have  approximately  a  simple 

(  1  &) 

quadratic  behavior  for  k  <  kp .  '  If  this  is  the  case  we 

find  from  equations  (2.28)  and  (2.29)  that  U  may  be 
evaluated  by  finding  the  diagonal  elements  of  G  at  an 
average  relative  momentum  k  =  /7"3  kp .  Then 
U  =  p<k|G|k>  and  the  average  energy  per  particle  is 


■  -=>  X  J 


33. 


E 


=  3  *1  k2 

5  2m  F 


+  ^p(kp)  <k | G I k> 


(2.30) 


A  similar  treatment  for  the  single  particle  potential 
is  possible.  Assume  that  we  may  replace  the  integration 
in  (2.24)  by  the  reaction  matrix  evaluated  at  an  average 
relative  momentum  defined  by  k2  =  ^(k2  +  0.6  kp  ).  Then 
U(kp)  is  approximately  quadratic  so  we  may  write  the  energy 
denominator  as; 

i 

e(k,k’  ,P)  *  k2  -k'2  +  2u[(P2  +k2f] 

*  k2  -k’2  +  2p(kp)  G(1^/P2Tk7+F7Fk^7  VP2+k2  +  0 .  bk2  ) 


Most  authors  treat  the  effect  of  the  center  of  mass 
momentum  by  replacing  P  by  its  average  over  the  Fermi  sea. 


=  kp(kF  -  k )  { 1  + 


k‘ 


} 


k  <  k. 


3kp(2kp+k) 


=  0 


k  >  kP. 


(2.31) 


We  have  used  the  value  P  =  v/7"3  kp  which  corresponds  to  the 
peak  in  the  distribution  of  the  number  of  pairs  of  particles 
with  center  of  mass  momentum  P.  Since  we  need  only  to 
evaluate  the  matrix  elements  <k|G|k>,  (with  k=  /TJ  kp) 
to  find  the  average  energy  per  particle,  the  energy 


denominator  is; 


. 

* 


34. 


e ( k , k ’ ,  P )  =  k2-k ' 2  +2 p ( kp ) <k | G | k> .  (2.32) 

(19) 

It  is  known  that  if  we  replace  the  energy 

denominator  by  a  free  particle  propagator  and  put  Q  =  1 
for  all  states,  then  the  diagonal  elements  of  the  reaction 
matrix  may  be  written  in  terms  of  the  phase  shifts. 

<k  |  G  |  k>  =  (2J+1)  A£ 

where  =  1/4  for  even  i  and  3/ 4  for  odd  l.  This  is  not 
a  particularly  good  approximation  for  nuclear  matter  as  it 
doesn't  account  for  the  observed  saturation^  ,  but  for 

i 

the  neutron  gas  at  low  densities  it  is  close  to  the  values 
of  the  reaction  matrix.  For  the  first  numerical  iteration 
of  the  problem  we  use  the  phase  shift  approximation  for  G 
in  the  denominator,  including  only  S-wave  contributions. 
Let  us  write  the  energy  denominator  as 


e ( k , k ' , P )  =  -{ y2 ( k)  +k'2}  (2.34) 


If  we  put  the  center  of  mass  momentum  identically  equal  to 
zero  and  neglect  the  exclusion  principle,  then  the  Green's 
function  may  be  evaluated  in  closed  form. 


G£(r ,r ' ) 


-2 


7T 


j  £(k'r)  j  £(k'r ' )  dk' 


o 


( k ’ 2+y2 ) 


r. . 


35. 


=  yj £( iyr  )  h£2 ^ ( iyr ’  ) 
=  yj  (iyr’ )  h^1 ) (iyr) 


for  r  <  r' 

(2.35) 

for  r  >  r  ’  . 


Then , 


U£(kr)  =  j  £ ( kr )  +y 


j£(iyd)  h( 1 ) (iyr)V(r ' )U£(kr ’ )r ’ 2 dr ' 


+  y 


j£(iyr)  h^^dyr’)  V(r')  U  (kr’)  r,2dr?. 


Put  U£  =  j  £ ( kr )  +  Ay  hj^dyr)  +  Byj£(iyr)  (2.36) 


The  boundary  conditions  are;  1.  U£(kc)  =  0 


2  .  U£d)  =  J  £  ( kr )  . 


These  imply  that  B  =  0,  and  A  =  ^  ^ 


yh£ 1  ( iyc ) 


;  or 


U£ ( kr )  =  j  £ ( kr )  +  hZ  ^iyr^  j£(kc). 


i 


(2.37) 


h£ 1 ^ ( iyc ) 


In  this  approxiamt ion ,  the  contribution  of  the  outer 
potential  to  the  reaction  matrix  is. 


<fJvIV 


{ j  £(kr )+j  £(kc ) 


c 


1  ^ ( iyr ) ' 

1 ^ ( iy° ) . 


}V(r )r 2dr. 


(2.38) 


•  *')  l  *  (“)  U  •  r:.  « 


This  is  identical  with  the  modified  Born  approximation  of 
reference  (20).  This  suggests  that  we  take  a  trial  wave 
function  to  be, 

U£T  =  J£(kr)-  h£  ^lx£r)  j£(kc)  (2.39) 

h[1 )(ix£c) 

where  x£  is  a  variational  parameter  that  will  depend  on 
angular  momentum  and  spin.  U£  has  the  form  of  a  plane  wave 
plus  a  term  corresponding  to  the  distortion  from  the  hard 
core . 

We  have  numerically  integrated  the  quadratic  functional 
J(U£T)  using  a  Gauss  Quadrature  integration  formula.  The 
contribution  of  the  outer  potential  to  the  diagonal  element 
of  the  reaction  matrix  was  determined  by  finding  the 
minimum  of  J  as  a  function  of  the  parameter  x£,  using  the 
Hamada-Johnst one  nuclear  potential  (see  page  45).  The 
first  iteration  was  started  with  the  phase  approximation 
for  the  reaction  matrix  in  the  energy  denominator.  The 
core  contribution  was  calculated  from  equation  (1.63). 

We  then  replaced  the  phase  shift  approximation  for  the 
reaction  matrix  in  the  denominator  by  the  calculated  G  and 
repeated  the  procedure.  This  process  converges  quite 
rapidly  for  the  density  range  .l<k0<2.7f  as  can  be  seen 

r 

from  the  graph  on  page  42.  The  integrations  have  been  done 
for  the  states  3S0,  3Po,  3Pi,  3?2.  The  total  energy  per 
particle  was  found  from  equation  (  2 . 30)  •  The  numerical 


r 


results  are  summarized  in  table  I  on  page  39  ,  with  the 

graph  of  energy  versus  density  on  page  40  .  We  note 

imediately  that  nuclear  forces  alone  are  not  sufficient  to 

produce  binding  of  the  neutron  gas  for  the  range  of  densitie 

( 23 ) 

considered.  M.  Barbi  '  has  reached  a  similar  conclusion 
for  the  three  neutron  problem  by  calculating  the  potential 
depths  needed  to  bind  three  neutrons  in  various  states 
using  a  numerical  solution  of  the  Euler-Lagrange  variational 
equations.  The  depths  obtained  are  not  compatible  with 
experiment  so  he  concludes  that  the  trineutron  doesn’t 
exist . 

Since  we  make  extensive  use  of  the  energy-density 
curve  later  on,  it  was  found  conveniant  to  express  the 
results  in  polynomial  form. 

E/N  =  a  j  kp  -  azkF  +  kF  (2.40) 

a3kz  +  a4 

Here  aj  is  just  the  kinetic  energy  term, 

3  W1 

ai  =i(7  IhT  =  12’41j9  a 2  =  .17763 


a  3 


a4  = 


10199 


.15972 


38. 


Several  other  authors  have  calculated  an  equation  of 

(  ^  rj  g  \ 

state  for  a  neutron  gas.  ^  5  55  '  Sood  and  Moszkowski  have 
treated  the  low  density  region  (k^  < . 5F~^ )  by  correcting 
the  free  pair  G  matrix  for  effects  of  the  Pauli  principle, 
assuming  a  separable  potential.  They  also  have  found  no 

( pi ) 

evidence  of  binding  by  nuclear  forces.  Salpeter  treats 

—1  —1 

the  moderate  density  range  ( . 6F  <  k„  <  1 . 9F  )  by 

r 

extrapolating  the  Weizsacker  semiempirical  mass  formula  for 

nuclear  matter  to  the  case  of  a  pure  neutron  gas.  His  values 

( 7  ) 

are  lower  than  ours.  Brueckner  et .  al.  have  used  the 

( 11 ) 

methods  developed  by  Brueckner  and  Gammel  to  evaluate 

self  consistant  single  particle  energies  and  total  energy 

using  a  Gammel-Thaler  potential.  Their  values  are  close 

to  ours  for  the  low  to  moderate  range  but  ours  fall  below 

theirs  for  higher  densities.  Levinger  and  Simmons v  '  have 

used  perturbation  theory  to  evaluate  the  energy  of  a  neutron 

gas  to  first  order  for  different  velocity  dependant  potentials. 

/  6  \ 

Sprung  et .  al .  '  have  used  the  Reid  potential  and  Brueckner- 

Bethe  many  body  theory  to  find  the  energy  of  a  neutron  gas 
and  have  also  estimated  the  fractional  density  of  protons 


in  the  system. 


» 
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Results  of  Numerical  Integration 
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Figure  1. 

Total  Energy  per  Particle  for  the  Neutron  Gas 
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45. 


The  Hamada-Johnston  Potential 


(1) 


The  model  consists  of  four  terms: 


V  =  V 


vts12 


vls(£ 


S  )  + 


vlll 1 2 


where  C,  T,  LS  and  LL  refer  to  central,  tensor,  linear  LS 
and  quadratic  LS  potentials  respectively,  and  L12  is  the 
operator  defined  by 


L  i  2 


( a  i 


a2)L: 


hi (ci -L) (a2 *L)  +  (ai -L) (a2 -L) } 


-*■  -> . 


=  {6lj  +  ( a i • a2 ) }L2  -  (L-S) 


The  VL(i  =  C,  T,  LS,  LL)  are  given  by, 

V  =  0. 08 (-y ) (ti *t2 ) (Si ’02 )Y(x)  { 1+a  Y(x ) +b  Y2 (x ) } 

C  j  O  O 

VT  =  0. 08(^-y  )  (ti  *t3  )Z(x)  (l+aTY(x)  +bTY2(x)} 

VLS  =  mGLS  Y2(x)  -t1+bLSY(x)> 

Vll  =  PGll  X”2  z(x)  {1+aLLY(x)  +bLLY 2  ^ x  ^ 

where  y  is  the  pion  mass,  and  x  is  measured  in  y  1 .  The 
numerical  parameters  are  on  the  next  page.  The  hard  core 
radius  is  x  =0.343  in  all  states  and  y=139*4  (Mev). 

Y ( x ) =  e'x  and  Z(x)=  (1+|  +  J  )Y(x) 


x 
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Chapter  IV 

The  Relativistic  Equations  of  Hydrostatic  Equilibrium 

To  date,  numerous  theoretical  models  of  the  origin 

and  evolution  of  stars  and  galaxies  have  been  proposed. 

These  may  be  considered  generally  as  falling  into  two  groups. 

In  one,  the  initial  stage  consists  of  diffuse  galactic 

matter  (mostly  hydrogen)  which  eventually  tends  to  condense 

with  the  lost  gravitational  energy  radiated  in  the  form  of 

light  and  heat.  The  temperature  of  the  interior  of  the 

mass  rises  until  nuclear  fusion  of  hydrogen  to  helium  occurs. 

The  next  period  in  the  life  of  the  star  is  the  so  called 

normal  period  during  which  all  the  stellar  thermonuclear 

processes  occur  and  is  also  the  longest  stage  in  the  life 

of  the  star.  Eventually  there  must  be  exhaustion  of  the 

nuclear  fuel  and  contraction  of  the  star  to  denser  states. 

What  happens  next  depends  specifically  on  the  configuration 

considered  and  will  be  discussed  later. 

In  the  second  group  of  evolutionary  models  the  initial 

state  is  a  superdense  body  with  the  normal  course  of 

evolution  proceding  to  states  of  lower  density.  It  is  not 

our  purpose  here  to  discuss  the  validity  of  these  models 

but  only  to  indicate  that  in  either  case  the  question  of 

matter  at  higher  than  normal  density  arises. 

It  has  been  suggested  that  a  system  of  baryons  at 

56 

the  endpoint  of  thermonuclear  evolution  form  Fe  in  a 

(24) 


body  centered  lattice. 


As  the  number  of  particles 
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increases,  the  gravitational  forces  increase  and  relativistic 

electrons  appear.  When  densities  exceed  about  3xl0^gm/cc 

all  of  the  heavy  ions  become  unstable  against  disintegration 

( 24 ) 

to  neutrons  by  means  of  electron  capture  and  the  system 

in  equilibrium  consists  mainly  of  neutrons  contaminated  with 

15 

some  protons  and  electrons.  Near  and  above  10  gm/cc  mesons 

and  other  baryons  appear.^0'1  For  the  range  lO^gm/cc  to 
15 

10  gm/cc  we  are  dealing  with  a  degenerate  gas  which  may  be 
considered  practically  as  consisting  only  of  neutrons. 
Cameron^^  has  considered  the  fraction  of  protons  and 
electrons  for  these  densities  and  has  shown  them  to  be 
about  three  orders  of  magnitude  lower  than  the  number  of 
neutrons . 

Many  works  have  predicted  the  existence  of  neutron 

stars.  Recent  discoveries  of  galactic  sources  of  X-rays 

renewed  hope  in  the  theories,  however  occultation  experiments 

have  shown  most  of  these  are  too  large  to  be  neutron  stars. 

( 25 ) 

As  early  as  1939,  Oppenheimer  and  co-workers  studied 

stars  during  the  last  stages  of  active  life  when  sources 
of  internal  energy  were  exhausted  and  the  star  had  cooled 
down.  They  predicted  equilibrium  configurations  of  neutron 
stars  at  zero  temperature  by  assuming  that  the  particles 
form  an  ideal  gas.  For  masses  greater  than  about  .7  solar 
mass  there  were  no  equilibrium  solutions.  When  thermonuclear 
sources  are  exhausted  heavier  stars  will  contract  unobstructed 
unless  they  can  loose  mass  by  some  mechanism.  It  is  now 
generally  accepted  that  supernovae  explosions  are  triggered 
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by  collapse  of  stellar  cores  followed  by  sudden  release  of 
gravitational  energy  or  ignition  of  nuclear  fuel  in  the 
remaining  envelope. 

More  recent  authors  have  studied  these  problems 
taking  into  account  in  various  ways  the  interparticle 
interactions  and  allowing  for  the  formation  of  hyperons  at 
higher  densities . } ^ ^ 3 ^  They  have  confirmed  the  result 
that  there  exists  an  upper  limit  around  one  solar  mass  for 
the  mass  of  a  neutron  star  beyond  which  there  is  no  stable 
configuration . 

The  behavior  of  nuclear  forces  at  extremely  high 

energies  is  known  only  approximately  and  the  interaction  of 

hyperons  is  not  known  at  all.  Consequently  it  is  difficult 

to  discuss  the  effect  that  formation  of  hyperons  has  on 

17 

superdense  configurations.  At  densities  around  10  gm/cc 

-1 2 

the  interparticle  distance  is  about  4x10  cm.  or  about 

the  radius  of  the  hard  core  for  nucleons.  Since  the 

-11 

interparticle  distance  is  of  the  order  of  10  cm  or  less, 

any  equations  should  include  the  structure  of  elementary 

( 29 ) 

particles  themselves.  Bahcall  and  Wolf  suggest  that  the 

concept  of  distinct  strongly  interacting  particles  is  not 
meaningful  for  densities  greater  than  about  eight  times 
nuclear  densities.  Also  the  distinction  between  Fermions 
and  Bosons  probably  disappears  around  this  density.  In 
this  case  the  star  should  be  considered  as  a  whole  and  not 
as  consisting  of  individual  particles.  It  seems  likely  now 
that  matter  should  follow  a  simple  asymtotic  equation  of 
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(31  30 ) 

state  at  higher  densities.  s  We  will  assume  this  type 
of  behavior  at  high  densities  and  attempt  to  describe  the 
low  to  moderate  density  range  correctly. 

In  the  case  of  out  own  sun,  the  ratio  of  the  gravitational 
radius  to  the  actual  radius  is  small  so  relativistic  effects 
are  small.  In  the  case  of  superdense  stellar  configurations 
however,  these  two  quantities  may  be  comparable.  In  other 
words  the  high  concentration  of  mass  produces  a  large 
curvature  of  space-time  near  the  star.  The  basis  of  our 
calculations  must  be  Einstein’s  gravitational  equations 
which  relate  the  energy-momentum  tensor  T  to  the 
fundamental  metric  tensor  g  .  In  the  system  of  units  for 
which  the  velocity  of  light  and  the  constant  of  gravitation 
are  equal  to  unity,  the  field  equations  are; 


-8ttT 


yv 


=  R 


l  -  ^R  g  +  Ag 
yi  v  toyv  ^yv 


(4.1) 


where  R  is  the  contracted  Riemann-Christof f el  tensor  and 
yv 

R  is  its  trace.  The  effects  of  the  last  term  are 

insignificant  within  the  size  of  the  solar  system  so  we 

may  consider  A  (the  cosmological  constant)  equal  to  zero. 

Physically  we  expect  a  static  spherical  distribution 

of  matter  so  we  choose  spatial  co-ordinates  r,  0,  <J>  and 

time  co-ordinate  t  to  reflect  this  symetry.  The  most 

general  static  line  element  exhibiting  spherical  symetry 

( 32 ) 

may  be  expressed  in  the  form  , 
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ds2  =  -ehr2  -  r2d02  -  r2sin20  dcf>2  +  evdt2  (4.2) 


where  A  and  v  are  functions  of  r  alone.  In  other  words 
the  components  of  the  metric  tensor  are; 


A  2 

gi  i  =  -e  g2  2  =  -rz 

g  3  3  =  -r  sin  0  g  4  4  =  e 

g  =  gyV  =  0  for  y  ?  v  and  gJ'J' 

M  v 


(4.3) 


( no  sum) . 


The  contracted  Riemann-Christof f el  tensor  R 
written  in  the  form; 


yv 


can  be 


R  =  {yv,a}  +  {ya  ,  $}{ v3 ,a}  +  — ^ - -log/-g  -{yv,a}-  log/-g. 

3xa  3xy3xv  3xa 

(4.4) 

where  g  =  |g  |  and  the  Christoff el  three-index  symbols 
with  curly  brackets  are  defined  by; 


{yv,a}  =  kg 


a  A 


9g  ,  ,  3g  ,  ,  9g 

toy A  +  tovA  +  toyv 


3x 


v 


3x 


v 


3x 


A 


(4.5) 


For  example; 


{11,1}  =  kg 


=  un 


3g i i  +  3g i i 


3g 


i  l 


3x 


3x 


3x  1 


i  dSl‘  =  =±.  1  (-eX)  =  A 

2gi  1  (  - )  A. A  3r 


3r 


2e' 


2 
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The  remaining  ones  are  got  in  a  similar  manner,  making  use 
of  the  relation  {yv,a}  =  {vp,a}.  The  non- vanishing  symbols 
corresponding  to  our  form  of  line  element  are;  (Equations 
95*5  of  reference  32) 


{11,1} 

=  h 

{31. 

,3} 

=  1/r 

{12,2} 

=  1/r 

{32. 

,3} 

=  ctn0 

{13,3} 

=  1/r 

{33: 

,1} 

=  -r  sin20e~^ 

{14,4} 

=  h  v' 

{33: 

» 2  } 

=  -sin0  cos0 

(4 

{21,2} 

=  1/r 

{41. 

,4} 

=  h  v' 

{22,1} 

-A 

=  -re 

{44, 

,1} 

1  v-A 
=  he  v 

{23,3} 

=  ctn0 . 

Using  these  and  equation  (4.4),  the  components  of  R  may 
be  calculated;  eg. 

Rn  =  ^|r{H,l>  +{11,1}2  +  {12,2}2  +  {13,3}2  +  {14,4}2 

+  —  log/-g  -  {11,1}  — log/^g  . 

9r 2  9r 

log/-g  =  \  log(e^+v  r4sin29) 

-^log/^g  =  4>(A"+v")  +  ~ 


, 
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—  log/— g  =  h(X"  +v")  -  “2 
3r 2  r 


R 


i  i  - 


-  ,,  "  2 


v  “  +v"  -X'v'  -X 

“T  ~  ~T~  - 


(4.7) 


Similarly , 


„  -A  ,  -A  f  v'-X'^  -! 

R 2  2  =  e  +  re  (  — - —  )  -  1 


(4.8) 


R33  =  rsin20  e  X  h(v'-X')  -  sin20  (1-e  X)  (4.9) 


R44  =  heV  X  i-v"  + 


v  '  A  ’  v  f  2  2v 


).  (4.10) 


2 


r 


It  will  prove  simpler  to  express  the  field  equations  in 
mixed  tensor  form; 

-87rTy  =  Ry  -  %R  gy  .  (4.11) 

v  v  °v 

Now  gy  =  6y  and  Ty  =  gyaT  .  Remembering  that  gyv  =  0 
Dv  v  v  av 

if  y  /  v,(4.11)  becomes, 


-8TrTy  =  gyy  R 
v  &  y  v 


sR 


6 


y  v . 


(4.12) 


The  non-vanishing  components  of  Ty  are  T} ,  T| ,  T3,  T4, 
which  may  be  found  by  substituting  the  previous  relations 


as  follows. 
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“8ttT  |  =  gnRir  3s(g11H11+  g22H22  +  g3  3R3  3  +  g 


or 


,  -16ttT|  =  — — 

1  gi  1 


1  _  Rll  -  R2  2  -  R  3  3  -  R  4  4 


gll  g  2  2  g  3  3  g  4  4 


-A  /  V  ’  2  .  v" 


-16ttT J  =  -e  A(  +  2 


A  '  v  '  _  AJ_  } 

r 


T 


+  lr 

r 


e  A+re  A^( v ’ -A ’ )  -  1 


1 


r2sin20 


rsin20  e  A^( 


-sin20(l-e  A) 


- 


-v 


e V  \(-2v"  +v  ’  A  ’  -  v’  2 


SttT;  =  -e’A{^-  +  ^2}  +  ^2 
r  r  r 


Similarly , 


8ttT2 


_-Arv"  -  A'v'  +  v'2  +  v’-A’t 
1 2  T~  ~T~  2r  1 


m3  _  m2 

3  x  2 


8^t; 


A ’  +1  N  ,  1 

e  (  —  -2  )  +  -2 

r  r  r 


Assuming  the  properties  of  a  perfect  fluid  (ie 
mechanical  medium  exerting  no  transverse  stresses), 
has  shown  (Reference  32,  equation  85.7)  that  we  can 
the  energy-momentum  tensor  in  the  general  form; 


Tyv  = 


( Po+Po ) 


dxy  dxv 
ds  ds 


g 


yv 


Po 


or , 


( Po+Po ) 


a  ,  v 
dx  dx 

ds  ds 


v 

gyPo 


4XJ 

u'-A’) 

\ 

)  . 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

,  a 

Tolman 

express 

(4.16) 


55. 


where  p0  and  p0  are  the  macroscopic  energy  density  and 
hydrostatic  pressure  respectively  measured  in  proper 
co-ordinates.  For  the  static  case  considered  here,  the 
components  of  fluid  velocity  are; 


dr  _  d6  d(j) 
ds  ds  ds 


0  — 

U  i  Ho 


dt  _  -v/2 

ds  e 


Putting  these  into  (4.16)  we  obtain; 


Ti  =  T2  =  T3  = 
12  3 


-Po  ,  Tj  =  (po+Po)  eve  -  p0  =  p0. 


v  -v 


We  now  have;  (droping  the  subscripts  on  p0and  p0) 


O  -A  r  v 1  +  1  -1  1 

ottp  =  e  —2  i  -  ~2 

^  2  r  r 


_  -ArV1  -  A  '  v  '  +  v'2  +  v ’-A ' , 

oirp  e  1 2  i(  Zj  2p  ^ 


(4.17) 


n  -  A  r  A  '  -1,1 

8ttP  =  e  {—  -2}+  -2 


Equate  the  first  two  equations  in  (4.17)  to  get; 


-A  rV" 

e  {  — 


-  A 1 v '  +  v'2  +  v ' -A ’  -  v’  -  1  t,  1  n  (4.18) 

! ~  ~TT  2r  ~  r2>+  r2  ~ 


Multiplying  be  2/r  and  rearranging,  this  is  equivalent  to; 


OJ  | 
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-X{A1  +  ^ 

r  r  2 


0  .  (4.19) 


Comparing  this  to  (4.17)  allows  us  to  write  (4.19)  as; 


(4.20) 


Einstein’s  field  equations  reduce  to  the  following  three 
equat ions ; 


8ttp 


(4.21) 


dp  =  -(p+p) v 
dr  2 


If  we  can  determine  A  and  v  as  functions  of  r,  the 

equilibrium  of  the  system  will  be  determined  by  adding 

some  equation  of  state  p  =  p(p).  In  the  empty  space 

surrounding  the  spherical  distribution  of  matter  p=p=0 

and  Schwarzschild ' s  exterior  solution  of  Einstein’s  field 

( 33 ) 

equations  is  valid. 


-A 


1  +  A/r 


e 


v 


B ( 1  +  A/r) 


e 


A  and  B  are  fixed  by  the  requirement  that  the  components 


57. 


of  the  metric  tensor  go  over  into  their  weak  field  form 
at  large  distances.  ie.  B  =  1,  A  =  -2m  where  m  is  the 
total  Newtonian  mass  of  the  matter  as  calculated  by  a 
distant  observer. 

The  boundary  of  the  sperical  distribution  is  at  r  =  rfe 
where  p  =  0  (remember  that  p>0  for  r<r^).  Using  the 
equation  of  state  e  =  e(p),  the  last  equation  in  (4.21) 
may  be  integrated. 


v(r)  =  v(rfe) 


p(r ) 

2dp 

.  P  +  p(p) 


or 


ev(r) 


y(  rh) 


exp 


p(r ) 

2dp 

, P+P(P) 

u 


(4.23) 


Now  ev  must  be  continuous  across  the  boundary.  From  (4.22) 


v(r)  2m  \ 

e  (1-  —  )  exp 

b 


r 

»  \ 

2dp 

1 

(P+P) 

J 

(4.24) 


e 


-A 


Introducing  a  new  variable  u(r)  =  3§r(l-e  ^)  or 
2  u. 

=  1  -  —  ,  the  second  equation  in  (4.21)  may  be  written. 


du 

dr 


4mp(r)  r2 


(4.25) 


The  first  equation  in  (4.21) 
e~^  =(1-  ~  )  and  using  v’  = 
equation . 


may  be  re-expressed  by  putting 

dp  (-2/(p+p))  from  the  third 
dr 
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dp  _  -(p+p)  {4-rrpr 3  +u}  (4.26) 

dr  r(r-2u) 


The  value  of  e'^r') 


must  be  continuous 


across  the  boundary. 


u(  r 


r,  n  -Mh)! 
=  _b  [  1-e  b  ] 

2 


=  ^b[l-(l- 
2 


M 


u(rfe)  is  the  observable  mass  of  the  star  as  determined  by 
its  gravitational  effect  on  a  distant  test  particle.  At 
this  point  we  make  a  change  in  notation  in  order  to  avoid 
confusion  later  on.  In  the  limit  of  large  density  we  have 
to  distinguish  between  the  total  energy  (mass)  density, 
(which  includes  the  sum  of  particle  masses  plus  energy 
of  motion  and  interaction,  not  including  gravitational 
energy)  and  the  density  of  the  rest  mass  of  the  particles. 
The  total  macroscopic  energy  density  previously  labled  p 
in  this  chapter  will  from  this  point  be  called  e.  The 
density  of  the  rest  mass  of  the  particles  will  be  p’=mp(k„) 

r 

where  p(kp)  is  the  number  density  of  particles  defined  in 
equation  (1.9) 5  and  m  is  the  rest  mass  of  the  neutron. 

The  system  will  then  be  described  by; 


du(r ) 
dr 


47rer  2 


dp  _  -(p+e)[  47Tpr3  +  u]  (4.27) 

dr  r(r-2u) 


e  =  e(p) 
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The  initial  conditions  will  have  the  form; 

u( 0 )  =0,  e  =  e(0)  =  e  ,  p  =  p(e(0))  (4.28) 

where  e  is  the  total  energy  density  at  the  center  of  the 
configuration.  Equations  (4.27)  are  integrated  simutaneously 
to  the  point  r  =  r^  where  p  =  0. 

The  proper  distance  in  a  gravitational  field  is 
determined  from 

drp  =  /_grr  dr  =  (1"  ^ ~±/ ?  (4.29) 

The  proper  mass  is  therefore 


dM 

_ P  = 

dr 


4  7T  p  ’  r 2 


r 

r-2u 

V.  / 


1/2 


(4.30) 


The  proper  mass  is  the  mass  the  star  would  have  if  the 
particles  were  dispersed  to  infinity.  The  binding  energy 
in  mass  units  is  the  proper  mass  minus  the  gravitational 
mass  and  is  given  by; 


dM 


b 


r 


1/2 

] 


dr 


47rr 2 


r-2u 


e 


(4.3]) 
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Chapter  V.  The  Maximum  Stable  Mass  of  a  Neutron  Star 


We  have  already  obtained  the  total  energy  per  particle 
of  the  neutron  gas  (E)  as  a  function  of  density  (or  k^)  in 

r 

equation  (2.40).  The  pressure  of  the  gas  is  given  by 

P  =  P2(|^)  =  1*80351  x  1031(ff  )  kp  dynes/cm2.  (5*1) 

P  Kp 

where  E  is  in  Mev  and  kp  in  F  ^ .  The  total  mass  density  e 
is  defined  by 


e  =  CpE(p)  +  p’c2]  =  [  E(Mev)  +  939*5]  kp  xO . 006011x10 1 3g/cc 

(5*2) 

where  kp  is  again  measured  in  F  ^ .  Remembering  that  we 
are  considering  a  static  structure  with  spherical  symetry 
and  also  that  p  and  e  are  in  proper  co-ordinates ,  the  stress 
energy  tensor  T^  may  be  written. 


T 


yv 


-p  0  0  0 

0  -p  0  0 

0  0  -p  0 

0  0  0  e 


(5*3) 


(30) 


Landau  and  Liftshitz 


state  that  for  a  non-interacting 


1  nl  _ji 

I 


6l . 


gas  the  trace 


of  T  must 
yv 


be  positive,  ie 


P  < 


(5.4) 


The  equality  sign  would  be  valid  only  for  a  photon  gas. 

( 11 ) 

However,  Zel’dovich  has  constructed  a  relat ivistically 

invariant  example  of  interaction  of  particles  with  a  field 
of  heavy  neutral  vector  mesons  which  yields  p  ■+  z  at  higher 
densities.  In  this  case  the  velocity  of  sound,  which  is 
given  by 


v 

s 


c/ 


de 


(5.5) 


is  always  less  than  or  equal  to  c,  the  velocity  of  light. 
This  seems  more  aesthetically  satisfying  than  the  previous 
example  which  yields 


v 

s 


(5.6) 


To  make  our  equation  of  state  have  the  correct  asymtotic 
behavior  at  high  densities  (p+e),  it  was  found  convenient 
to  modify  equation  (2.40)  as 


E  = 


v 


a  ik 


2 

F 


a2kp  +  kp 
a3kp  +  a 4 


X(kp)  +  [  1  -  x(kF)]k|  (5.7) 


'* 
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where 


X(kp)  =  (  1  +  exp[  kj,  -  8  ]  ) 


(5-8) 


This  leaves  the  low  to  moderate  density  regions  unchanged 
but  at  £  ^  1016g/cc,  p  e  .  Ours  then  is  a  parametric 
equation  of  state. 


e  =  e(kp) 
P  =  p(kp) 


(5.9) 


which  may  be  compared  in  Figure  6  with  those  of  Cameron ^ ^ ^ . 

We  have  used  this  equation  of  state  to  integrate  the 

equations  of  equilibrium  (4.27),  (4.30),  and  (4.31)  for 

various  central  energy  densities.  The  results  are 

summarized  in  table  III  and  Figure  7. 

The  quantities  in  equations  (4.27),  (4.30),  and  (4.31) 

are  in  relativistic  units  (c  =  G  =  1  and  M4  c  5/32tt26i  3  =  K 

=  1/4tt).  In  order  to  change  back  to  c.g.s.  units  we 

(  34 ) 

must  multiply  by  the  following  conversion  factors. 


3 

length:  r0  =  2/2?(tf/Mc)2  (c//GM) 

=  13.7  km  -  2  x  10_5R  (5.10) 

© 

mass:  m0  =  r0c2/G  =  1.85  x  1034g  =  9-29  M 

G) 

pressure:  po  =  (M4  c  5  /  3  2  tt  2tf  3 )  4tt  =  6.46  x  10  3  6dynes/cm2 
density:  p0  =  Po/c2  =  7-15  x  1015g/cm3. 
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The  neutron  gas  is  relativistic  in  the  core  of  some 
of  the  configurations  considered.  Strictly  speaking  we 
should  have  used  the  relativistic  expression  for  the  kinetic 
energy , 


K.E. 


c  (m  2c  2+p 


2)  1/2 


(5.11) 


instead  of 

K.E.  =  (me 2+  p2/2m)  (5-12) 

(  35 ) 

Saakyan  has  integrated  the  equations  of  equilibrium 

using  equations  of  state  for  real  neutron  gases  with  either 
(5-11)  or  (5-12)  and  has  shown  that  the  mass  radius  values 
agree  within  three  significant  figures.  The  reason  for  this 
is  that  the  repulsive  forces  become  dominant  at  energies 
where  the  neutron  gas  is  relativistic  so  it  makes  little 
difference  which  of  (5*11)  of  (5-12)  is  chosen. 

The  first  qualitative  feature  that  appears  when  we 

plot  the  observable  mass  against  the  central  energy  density 

(Figure  7)  is  the  existence  of  a  maximum  in  the  curve 

at  £  =  6 . 3  x  10  g/cm  and  M  =  1.6Mn.  This  mass  maximum 

c  w 

is  commonly  called  the  OV  (Oppenheimer  and  Volkoff)  maximum. 
As  previously  mentioned,  this  maximum  was  first  noticed  by 
Oppenheimer  and  Volkoff  in  1939  when  they  integrated  the 
equilibrium  equations  using  the  equation  of  state  for  an 
ideal  Fermi  gas . 
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The  appearance  of  a  maximum  in  the  curve  is  not 
unexpected  since  it  is  easy  to  show  that  the  mass  of  an 
equilibrium  configuration  cannot  be  arbitrarily  large. 

It  is  well  known  that  the  gravitational  force  on  a  test 
mass  goes  to  infinity  if 

2GM 

r  r  =  -  .  _  _  _ . 

g  c  2  (5.13) 


Here  r^  is  called  the  gravitational  radius,  G  is  the 
gravitational  constant  and  M  is  the  observable  mass  of  the 
configuration  considered.  Define  an  average  energy  density 


(5.14) 


then  R  >  r  implies 
g 


M  < 


(2g)3/2(|  >1/2^72 


(5.15) 


The  largest  possible  mass  of  the  equilibrium  configuration 
must  decrease  with  increasing  e. 

(of.) 

Zel ' dovich v J  ;  has  done  an  analysis  of  the  exact 
equilibrium  conditions  for  the  star  and  has  shown  that  a 
configuration  is  stable  if 


dM  >  0.  (5.16) 

de 

c 


This  is  a  natural  criterion  since  in  the  stable  state,  the 


' 


. 
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the  addition  of  mass  causes  compression  with  an  increase 
in  pressure  to  compensate  the  gravitational  force. 

Neutron  stars  whose  observable  mass  is  less  than  1.6  M 

© 

will  be  dynamically  stable.  The  star  as  a  whole  will  not 

consist  only  of  neutrons  since  in  the  outer  regions  the 

pressure  is  not  sufficient  for  the  existance  of  stable 

neutrons.  It  was  apparent  during  the  integrations  that  the 

density  in  the  interior  of  the  star  remained  remarkably 

constant  for  the  greatest  part  of  the  star  radius  and  then 

fell  rapidly  to  zero.  This  is  the  same  observation  that 
( -34 ) 

Cameron  has  made.  He  has  shown  that  models  of  low 

density  may  have  extensive  electronic-nuclear  envelopes 

but  for  higher  central  energy  densities  (  1014g/cm3)  the 

shells  diminish  and  are  no  longer  important.  Cameron  also 

has  integrated  the  equilibrium  equations  for  pure  neutron 

configurations  as  well  as  for  models  in  which  the  composition 

changes  with  density.  The  effect  of  changes  in  composition 

on  the  mass  radius  values  is  very  small  compared  to  the 

effect  that  different  assumptions  about  nuclear  forces  has. 

We  feel  completely  justified  in  regarding  the  star  (for  the 

range  10 14  <  e  <  1016g/cm  )  as  being  composed  of  a  neutron 

0 

core  surrounded  by  thin  or  negligible  envelopes  of  electrons 

(  ”5  7  ) 

and  nuclei.  Chiu  and  Salpeter  have  determined  the 

surface  temperature  and  luminosity  for  neutron  stars.  The 
surface  temperature  will  fall  in  the  range  106  ■+  107  °K. 


For  temperatures  even  as  high  as  109°K  and  e  =  1015g/cm3 
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the  ratio  of  Fermi  energy  to  thermal  energy  is  about  900 

so  the  cold  matter  approximation  is  justified.  For 

temperatures  any  higher  than  this  neutrino  processes 

would  rapidly  dissipate  energy.  Also  the  high  thermal 

conductivity  of  the  degenerate  neutron  gas  and  small  dimensions 

of  the  neutron  star  insure  that  the  temperature  is  low. 

The  luminosity  is  thousands  of  times  greater  than  our  sun 

but  the  maximum  of  the  energy  distribution  lies  in  the  X-ray 

region.  Most  of  the  radiation  will  not  penetrate  our 

atmosphere  so  discovery  of  a  neutron  star  may  have  to  wait 

(9) 

for  an  orbiting  telescope.  Chiu  has  said  that  in  order 
to  be  visible  on  the  surface  of  the  earth  by  the  largest 
telescope  today,  a  neutron  star  would  have  to  be  within 
our  own  solar  system.  It  is  no  small  wonder  that  a  neutron 
star  has  not  been  discovered  to  this  date. 

At  this  point  we  must  say  something  about  the  question 
regarding  the  behavior  of  configurations  with  masses  greater 
than  the  OV  maximum  for  which  there  are  no  stable  solutions 
to  the  relativistic  equations .  To  either  side  of  the  OV 
maximum  the  gravitational  force  exceeds  the  pressure  and 
causes  the  star  to  contract.  To  the  low  density  side 
however,  the  pressure  increases  and  returns  the  star  to 
equilibrium  while  to  the  high  density  side  the  forces  take 
the  star  further  from  equilibrium  and  the  star  contracts 
more.  When  an  amount  of  cold  matter  is  assembled  with  a 
mass  exceeding  the  critical  figure  of  1.6  M  ,  then  unless 

Q 
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(38) 


it  can  shed  mass  by  some  mechanism  it  will  become  unstable 
against  collapse  and  contract  unobstructed.  This  is  the 
phase  known  as  gravitational  collapse.  It  is  beyond  the 
scope  of  the  present  work  to  attempt  any  calculations  on 
these  processes  but  for  the  sake  of  completness  we  will 
briefly  mention  some  of  the  current  ideas.  M.  A.  Podurets 
has  shown  that  the  pressure  forces  are  not  important  for 
calculating  the  dynamics  of  collapse.  For  a  first 
approximation  it  is  possible  to  consider  the  free  fall  of  a 

(  QQ  ) 

particle  in  a  Schwart zchild  field.  Zel'dovich  and  Novikov 
explain  the  process  as  follows.  A  remote  observer  sees  the 
star  collapsing  at  a  large  rate  when  its  radius  is  still 
much  larger  than  r  .  The  velocity  of  collapse  goes  to  zero 

o 

as  r  -*  r  (as  seen  by  the  remote  observer).  The  time  taken 
g 

to  fall  to  r  as  seen  by  a  local  observer  is  finite  however 
g 

and  the  star  continues  to  collapse  after  reaching  r  .  At 
r  =  r^  the  Schwarzchild  singularity  occurs.  On  this  surface 
the  curvature  of  space  becomes  infinite  (no  light  emitted 
from  this  surface  will  reach  us).  As  seen  by  the  external 
viewer,  when  the  radius  reaches  r^  the  star  stops  radiating 
almost  instantaneously  (although  the  star  still  interacts 
with  surrounding  bodies  through  its  own  gravitational  field). 
The  star  has  effectively  cut  itself  off  from  the  rest  of  the 
universe  by  producing  such  a  large  curvature  that  not  even 
photons  may  leave  it . 
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At  the  International  Symposium  on  Gravitational 
Collapse  and  Relativistic  Astrophysics  in  Dallas  ( 1 9 6 3 ) 
it  was  concluded  that  there  was  no  principle  which  could 
prevent  the  collapse  from  proceeding  to  completion  and  in 
the  process  the  constituent  baryons  would  disappear. 

This  is  a  violation  of  present  day  ideas  and  may  represent 
the  appearence  of  new  physical  principles  or  may  in  fact 
indicate  that  we  must  find  some  other  approach  to  describe 
the  collapse.  At  any  rate  the  problem  as  it  now  stands 
is  unresolved. 
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Table  III.  Characteristics  of  Neutron  Star  Models. 
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Figure  7 .  The  Maximum  Stable  Mass  of  a  Neutron  Star 
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